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On the basis of relativistic mean field calculations, we predict that the spin-orbit splitting of p 3 / 2 
and P1/2 neutron orbits depends sensitively on the magnitude of the proton density near the center 
of the nucleus, and in particular on the occupation of s 1 / 2 proton orbits. We focus on two exotic 
nuclei, 46 Ar and 206 Hg, in which the presence of a pair of Sj/2 proton holes is predicted to cause the 
splitting between the p 3 / 2 and Px/ 2 neutron orbits near the Fermi surface to be much smaller than in 
the nearby doubly-magic nuclei 48 Ca and 208 Pb. We note that these two exotic nuclei depart from 
the long-standing paradigm of a central potential proportional to the ground state baryon density 
and a spin-orbit potential proportional to the derivative of the central potential. 
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One of the primary motivations for the study of ex- 
otic nuclei is to search for novel shell structure ef- 
fects. A large amount of attention has been paid to 
the possibility that the spin-orbit force on high-j neu- 
tron orbits weakens in nuclei near the neutron drip line 
HaSISiSIIQSmEEQ The neutron 
magic numbers for stable nuclei rely on the effect of the 
strong spin-orbit force on high-j orbits, so the weaken- 
ing of this force has the potential to change the neutron 
magic numbers in neutron-rich nuclei. The possibility 
of the narrowing or collapse of the N = 28 major shell 
closure in neutron-rich nuclei near Si has attracted con- 
siderable attention because these isotopes are becoming 
accessible to experiments 0, ^| . The two most 
important reasons generally given for the decline of the 
spin-orbit force on high-j neutron orbits near the neutron 
drip line are the large neutron surface diffuseness and the 
influence of the continuum in these nuclei (T^J |2(j ■ 

In the present communication, we predict a novel shell 
structure effect having to do with spin-orbit splitting in 
low-j neutron orbits — namely, p orbits. The dramatic 
decrease in the spin-orbit splitting described here is not 
caused by the neutron density near the nuclear surface, 
but rather by the proton density in the nuclear interior. 
The two nuclei for which we make specific predictions, 
46 Ar and 206 Hg, are exotic but within two protons of the 
valley of stability. We make these predictions using the 
relativistic mean field theory, which has also been used 
to study the s pin- orbit splitting of high-j orbits in exotic 
nuclei 0, IE H1HG1 

The relativistic mean field calculation reported here is 
identical to the calculation used in Ref. [2l| to predict the 
properties of neutron-rich nuclei over a wide mass range. 
The model used in Ref. |21l is based on a Lagrangian 
developed in Refs. p2l l23j that includes novel nonlin- 
ear couplings between the isoscalar and isovector mesons. 
These new terms, which supplement the phenomenologi- 
cally successful Lagrangians of Refs. |2J, |2Jj, |2f| , modify 
the density dependence of the symmetry energy with- 
out changing ground state properties that are well es- 
tablished experimentally. Modifications to the poorly 



known density dependence of the symmetry energy in- 
duces interesting correlations between the neutron skin 
of heav y n uclei and a variety of neutron-star proper- 
ties miiiiiii. 

In both doubly-magic nuclei 48 Ca and 208 Pb, the high- 
est lying proton orbits below the Fermi surface (or the 
lowest energy proton hole states in 47 K and 207 T1) are 
Si/2 orbits. The effect of removing a pair of S1/2 pro- 
tons from 48 Ca and 208 Pb is illustrated in Fig. ^ which 
compares the proton densities of 46 Ar and 48 Ca (up- 
per panel), and the proton densities of 206 Hg and 208 Pb 
(lower panel). It is important to note that the predicted 
root-mean-square charge radii for 48 Ca and 208 Pb are 
in excellent agreement with experiment |29| . As the Si/ 2 
wavefunctions are strongly peaked in the center of the nu- 
cleus, the removal of these protons from 48 Ca and 208 Pb 
results in sharply reduced proton densities in the centers 
of 46 Ar and 206 Hg. This, in turn, causes a sharp increase 
in the magnitude of the spin-orbit interaction in the nu- 
clear interior. Figure Ufa) illustrates this effect in 208 Ph 
and 206 Hg; inside of 2 fm, V so is much stronger — and 
of the opposite sign — in the Sy 2 nucleus 206 Hg than in 

the doubly-magic 208 Pb core. This unconventional be- 
havior of the spin-orbit potential is intimately related to 
the Lorentz structure of the Dirac mean fields. While the 
depletion of S\/2 proton strength manifests itself in both 
the vector and scalar densities, the "proton-hole" disap- 
pears from the central potential as a result of the sensitive 
cancellation between the attractive scalar and the repul- 
sive vector potentials [30|. In contrast, (the derivatives 
of) the scalar and vector potentials add constructively in 
the spin-orbit potential and the development of a non- 
trivial spin-orbit structure in the interior of the nucleus 
ensues. 

Figures|3b)-(c) illustrate the effect of folding the spin- 
orbit potential in 208 Pb and 206 Hg with the 2>P\/2 and 
3p 3 / 2 neutron wavefunctions. That is, we display 
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FIG. 1: Proton (point) densities for (a) 46 Ar and 48 Ca 
and for (b) 206 Hg and 208 Pb computed using the relativistic 
parametrization of Ref. |24^ . The development of a "proton 
hole" in the interior of the nucleus is readily observed. 
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FIG. 2: (a) The Schrodinger-equivalent spin-orbit potential 
for 208 Pb (solid line) and 206 Hg (dashed line). Panels (b) and 
(c) display the effect of folding the spin-orbit potential with 
the Schrodinger-equivalent p-orbitals, as defined in Eq. Q. 
The arrows point to a first-order estimate of the spin-orbit 
splitting. 



Expressions for the Schrodinger-equivalent spin-orbit 
potential and wavefunctions [u{r)\ may be found in 
Ref. [33. For the purposes of this study, normalized 
wavefunctions have been used. Note that the above quan- 
tity, while not exact, provides an accurate (first-order) 
estimate of the P3/2~Pi/2 spin-orbit splitting AV so = 
AV so (r->oo). 

The combined effect of a strong increase in V so in the 
interior of 206 Hg together with neutron wavefunctions 
that are much larger at small radii than larger j orbits, 
yields a big effect on the integrated spin-orbit energy and, 
therefore, on the spin-orbit splitting for the p-neutrons. 
Indeed, this effect leads to the collapse of the P3/2-P1/2 
spin-orbit splitting: from -0.60 MeV in 208 Pb to -0.02 
MeV in 206 Hg. A similar effect occurs in 46 Ar relative 
to 48 Ca: the spin-orbit splitting is reduced by almost an 
order of magnitude, as shown in Fig. [3] Figure 0] displays 
the exact values. We conclude that in the two exotic nu- 
clei 46 Ar and 206 Hg, the spin-orbit interaction ceases to 
be a surface-dominated phenomenon. 

Figure0]summarizes the calculations and compares the 
experimental and calculated binding energies for £3/2,1/2 
neutron orbits in 48 Ca and 208 Pb. It should be empha- 
sized that while models with nonlinear couplings between 



the isoscalar and isovector mesons change the energy of 
the individual p 3 / 2 and P1/2 orbitals slightly, the predic- 
tion for their spin-orbit splitting is largely model inde- 
pendent. Thus, we limit ourselves to the original NL3 set 
of Ref. 0]. The experimental binding energies for 48 Ca 
are taken from the 48 Ca(c?,p) measurement of Uozumi et 
al. l32| and the mass compilation of Audi and collabora- 
tors [33. The 208 Pb(p,d) data used to extract the bind- 
ing energies for 208 Pb are taken from the compilation of 
Martin [34|]; the mass data are taken from Ref. |33[ . 

For stable nuclei, the standard experimental technique 
for mapping single neutron strength in a nucleus is to 
use a stripping reaction such as (d,p). To differentiate 
between spin-orbit partners (such as P3/2 and Pi /? ) a po- 
larized deuteron beam would be used (as in 32]). For 
the exotic nucleus 46 Ar, such a measurement would be 
performed in inverse kinematics with a 46 Ar beam and 
polarized deuteron target. The measurement would fur- 
ther be complicated by the likelihood that the P3/2 and 
P1/2 strengths would be somewhat fragmented, as they 
are in 51 Ti, 53 Cr and 55 Fe |H. In 206 Hg, the £3/2,1/2 or- 
bits would be observed as holes, requiring the use of the 
pickup reaction (p,d) in inverse kinematics, once again 
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FIG. 3: (a) The Schrodinger-equivalent spin-orbit potential 
for 48 Ca (solid line) and 46 Ar (dashed line). Panels (b) and 
(c) display the effect of folding the spin-orbit potential with 
the Schrodinger-equivalent p-orbitals, as defined in Eq. fl). 
The arrows point to a first-order estimate of the spin-orbit 
splitting. 



with a polarized target to differentiate between spin-orbit 
partners. For example, the normal kinematics experi- 
ment 208 Pb(p,d) with a polarized beam is reported in 

it 

In summary, we have used relativistic mean field cal- 
culations to predict that the spin-orbit splitting of p 3 / 2 
and pi /2 neutron orbits depends sensitively on the magni- 
tude of the proton density near the center of the nucleus, 
and in particular on the occupation of s±/2 proton orbits. 
The quenching (or collapse) of the spin-orbit splitting in 
high-j neutron orbits has been advertised as the hallmark 
for novel nuclear-structure effects in neutron-rich nuclei. 
This collapse is associated with the development of a dif- 
fuse neutron-rich surface. In this communication we have 
proposed a new mechanism for the collapse of the spin- 
orbit splitting — but among low-j neutron orbits. This 
mechanism is based, not on a rearrangement of the neu- 
tron density at the surface of the nucleus, but rather, on 
a depletion of the proton density in the nuclear interior. 
Two exotic nuclei, 46 Ar and 206 Hg, may be accessible 
for the study of this effect. In these nuclei we predict 
that the presence of a pair of Sim proton holes causes 



the splitting between the P3/2 and P1/2 neutron orbits 
near the Fermi surface to be much smaller than in the 
nearby doubly-magic nuclei 48 Ca and 208 Pb. Thus these 
two exotic nuclei, only two protons away from being dou- 
bly magic, deviate from a long-standing paradigm that 
has been applied with enormous success in both struc- 
ture and reaction calculations, namely, that of a central 
potential proportional to the ground state baryon density 
and a spin-orbit potential proportional to the derivative 
of the central potential. Finally, while the effects pro- 
posed here are likely to be modified by correlations that 
go beyond the mean-field approximation, we trust that 
most of its novel qualitative features will hold true. 
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FIG. 4: Comparison between the experimental and calculated 
P3/2-P1/2 spin-orbit splitting for the doubly-magic nuclei 48 Ca 
and 208 Pb. Also shown is the predicted collapse of the spin- 
orbit splitting in the two exotic nuclei 46 Ar and 206 Hg. 
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